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Abstract. In recent years light-front quantisation has been extended to allow for 
a consistent treatment of systems at finite temperature and density. This is in 
particular interesting for an investigation of the processes in nuclear matter under 
extreme condition as occurring, e.g., during a heavy ion collision. Utilising a Dyson 
expansion to the TV-point Green functions at finite temperature and density we focus 
on the occurrence of pionic and scalar diquark dynamics in quark matter and compute 
the masses and the Mott dissociation using a separable i-matrix approach. For 
the scalar quark-quark correlation we determine the critical temperature of colour 
superconductivity using the Thouless criterion. On the same footing the properties of 
the nucleon in a medium of quark matter are computed within a Faddeev approach. 
Critical lines for nucleon breakup are given. Presently, we use a light-front Nambu- 
Jona-Lasinio model that allows us to compare these results of this novel approach to 
the more traditional instant form approach, where applicable. 
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1. Introduction 

Light-front quantisation recognized by Dirac [1] provides a framework to describe the 
perturbative and the nonperturbative regime of quantum chromodynamics (QCD), for 
an overview see e.g. j2]. The importance of light-front quantisation as a complement to 
Monte Carlo simulations has also recently been emphasized by Ken Wilson [3]. 

A description involving both regimes is necessary, if one is interested in the dynamics 
close to the chiral and the confinement-deconfinement transition that is suggested 
by lattice QCD gl U U E], model analyses [9j M, HI E2], and through recent 
interpretations of certain RHIC results [13l dH [151 E3 E3 EEl EH] . It seems that close 
to phase transition the quark gluon plasma does not appear as a weakly interacting 
gas of quarks and gluons but rather as a strongly correlated system. In the framework 
of many-body Green functions [20l [21], the inclusion of correlations can be achieved 
by the Dyson equation approach developed in the context of nonrelativistic nuclear 
physics, reviewed in [22]. This approach has been generalised to light-front quantisation 
to comply with the requirements of special relativity in |23j, where also the necessary 
ingredients such as distribution functions, Green functions and Matsubara frequencies 
have been given. An early exploration of finite temperature QED(1+1) in discrete light 
cone quantisation has been done in Ref. [23] and the same theory has been reexamined 
in [25]. Further developments have been carried out and more examples given in Refs. 

[261 [271 [281 [291 [301 [311 [321 [Ml [Ml 

A formal framework of covariant calculations at finite temperatures in instant form 
has been given in Refs. [36l|37ll38] in a different context. Combination of the covariant 
statistical physics and light-front quantisation leads to the notion of general light-front 
coordinates first introduced in |30j. These coordinates allow the formulation of thermal 
field theory within the light-front frame and help to overcome the feature that the 
heat bath is moving with time-like velocity, while the time direction vector is light-like. 
Choosing a certain set of parameters the oblique frame in [30] reduces to the general 
light-front frame suggested in [23j |27] which correspond to a heat bath at rest in the 
instant form. The transformation is repeated here for completeness [30J 

X° = X + , Xq= x , x ± = (x 1 ,x 2 ), (1) 

,-y*3 ,-y»3 ,->-» — Q^yi nr*— L ( ry 1 rf* ^ A 

Ju JU y JU *J £iJU - JU 1 JU - tAJ j . 

One recognizes that in (EJ) light cone and instant form coordinates get mixed and 
the covariant time component is light cone time while the contravariant component 
is ordinary time. The metric can be derived by requiring invariance of the Minkowski 
line element and the important condition goo > holds. Therefore the medium velocity 
can be expressed as = (1/^/ggg, 0, 0, 0) and fulfills u 2 = 1. Only few applications of 
thermal field theory in the general light cone frame are found in the literature. 

Here we explore finite temperature and density properties of few-body, namely one-, 
two- and three-body systems in the Nambu-Jona-Lasinio (NJL) model that reproduces 
some low energy phenomenology of QCD, in particular its chiral properties, in a rather 
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transparent way, and yet includes spontaneous symmetry breaking as a nontrivial feature 
to be handled by light-front quantisation. The investigation of the one- and two- 
body sector confirms known results, like the phase boundary of chiral restoration, the 
pion Mott transition, the onset of two flavour colour superconductivity, from earlier 
calculations [39J HO] and exemplifies the equivalence of thermal field theory in the front 
and instant form. The three-body part consists of new results within this approach. 

This paper is organized as follows. In the second Section we derive the explicit 
form of the meanfield thermal fermion propagator using general light cone coordinates 
and thereby showing that one recovers the propagator expected in the ordinary light- 
front frame. In Section 3 the propagator is embedded into the general framework of the 
Dyson expansion of many-body Green functions. A systematic treatment for two- and 
three-body correlations in the quark matter medium is developed. The last two Sections 
are concerned with applications of the framework to correlations of two and three light 
quarks in meanfield approximation. The Mott dissociation lines of the mesons and 
nucleons are computed and plotted in the phase diagram of quark matter. 



2. Light-front equilibrium Green functions 

Let us summarise aspects of thermal field theory in the light-front form and derive the 
free, in-medium Green function. In the generalised coordinates (pQ) the product u ■ P is 
given by Pg/ \fg^, and the grand canonical partition function reads 

Z G = Tr exp {-(Po/Vm " t*N)/T} = Tr exp {-(P + + P_ - fiN)/T} , (2) 

where T is the scalar rest frame temperature and \i the chemical potential belonging 
to the charge N. We introduce H g q = P + + P_ — fiN and the statistical operator 
pz = e~ HcS ^ T /Zq. The four-momentum operator and the charge are 

P»= J d 2 x ± dx* T^(x) = / d 2 x ± dx 3 T+"(x), (3) 
N = [ d 2 x ± dx* f°{x) = [ d 2 x ± dx 3 J + {x), (4) 



where T^ v (x) denotes the energy momentum tensor and J^(x) the conserved current. 
Note that in Q and (jl]) the light cone densities are integrated over three dimensional 
space and not the light-front plane. 

The general light cone time-ordered (causal) Green function is 

G a0 (x -y) = 6(x Q - y«) G> p (x - y) + 6(y° - x°) Q< p (x - y), (5) 

with the correlation functions 

g>^-y\ X -y)= -i(y a (x)Vp(y)), (6) 
Q<^-y\ X -y)= T (-i)(Vp(y)y a (x)), (7) 

and the notation 3 ). Here and in the following the upper (lower) sign is for 

fermions (bosons). The anti-causal Green function can be defined correspondingly. In 
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the generalised Heisenberg picture the light cone time dependence of operators is given 
by 

0(x\ x) = e iHrffX °£>(0, x)e- iH ° sX \ (8) 

This differs from the regular Heisenberg picture of the vacuum by the thermodynamic 
constraints. In equilibrium the average is taken over the (equilibrium) grand canonical 
statistical operator pc, viz. (•■■) = tr{pc ■ ■ ■} and in addition and Q>a are related 
by (anti) periodic boundary conditions, 

g< (x°,x) = ±g> (x°-zp,x). (9) 

Therefore, in equilibrium, only one Green function needs to be considered, and, 
alternatively to (j3j), we may introduce the thermodynamic Green function [21] 

07/ = -<Tr*«(r)^(70>. (10) 

Setting x + = —ir ( "imaginary time" ) in (jSJ) we achieve the thermodynamic Heisenberg 
picture. The Green functions (JSJ) and fflUj) are related by their spectral representation 
utilising analytic continuation of the respective spectral functions just as in the instant 
form case [21]. To be more specific, in momentum space we introduce the spectral 
function A(ko,k) 

G < (k- ,k) = f(k- ,k)A(k- ,k) (11) 
Q > (h,k) = (l~f(h,k))A(k- ,k) (12) 

where A(ko,k) = i(G > (kQ,k) — G < (kQ, k)) and k abbreviates (hi, k^, fcj). Using <^ along 
with (JTTI) and ( Tl2l) the causal Green function can be represented as 



dko 
2(2tt) 



f(k- ,k)A(k- ,k) | (l-/(fc 5 ,fc))A(fc g ,fc) 

2^0 2 ^ 2^5 2 ^ 



(13) 



A similar equation holds for (flOl) . where the is term can be dropped 



g(kik)= f (14) 

y{ J 2(2tt) l]fc» - ifo 1 J 

— oo 

and the Matsubara frequencies are given by 

{i(2n + l)irT + a fermions, 
(15) 
i2nirT + p bosons. 

One may introduce traditional light cone coordinates via (pQ) into the equations ffTTj) to 
(TT4j) . since they only depend on the difference fcg — Then the Matsubara frequencies 
in A; - read 

{i(2n + lWT + u — 1/e + fermions, 
. (16) 
osons, 
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and have been given before in Ref. [23] without the introduction of the general light- 
front frame. For an ideal gas (and also in Hartree-Fock approximation utilised later on) 
the spectral function is 

A(kQ, k) = 2ix (jk + m) e(ko)5(k 2 - m 2 ) 

= 2n^^e(k- 3 )S(h-k- , on ), (17) 

where e(x) denotes the sign function. Inserting (TlTj) into ({TBI and transforming from 
general light cone frame back to the traditional one leads to 

rn , lKn + m f f(k + ,k ± ) 1 -/(£;+, fcj 



k ^on ^on 

with k~ n = (k\ + m?)/k + . Separating positive from negative k + components and 
introducing the grand canonical Fermi distribution functions of particles f + = f and 
antiparticles f~(k + ) = 1 — f + (—k + ) 



Eq. ffl8l) changes to 



(19) 



~ k~ — k~ n + is k+ (W ^ + jfe-_jfe-_ ie fc+ ' (2U) 



+ 7fc on + m 6(-k + ) ^_ + 7fc on + m 9(-k + ) _ 

which has been given before in Ref. [UJ following naively a direct approach to light-front 
Green functions. However, the results are as shown equal and the particle propagator, 
that is setting /" = in (TSUI) , can be found in Ref. [23]. In comparison to (1201) the Green 
function given in Ref. [271 [32] were derived using the closed time path formalism and have 
therefore 2x2 matrix structure. These Green function are suited to non-equilibrium 
situations, while we here concentrate on the equilibrium or close-to-equilibrium systems, 
which is why ( 12U1) is sufficient. 



3. Dyson expansion 

To investigate bound states in hot and dense quark matter we use techniques of the 
many-body Green functions organizing the equations into a Dyson expansion that 
leads to a hierarchy of linked cluster equations. This approach for the investigation 
of correlations in many-body systems at finite temperature and density was derived by 
P. Schuck and collaborators [22] . Furthermore this formalism was extended to the light- 
front in [23] to investigate three-quark correlations in quark matter. In the following 
we present a covariant derivation of the Dyson approach on the light-front and its 
application for mesonic and baryonic bound states. The Dyson approach to many- 
body Green functions in the light-front quantisation allows to calculate systematically 
the properties of few-body clusters, in particular of the two-quark bound states (viz. 
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pion) and of the three-quark clusters (viz. nucleon). This sets up the framework to 
investigate the change from nuclear to quark matter. One starts from the generalisation 
of the casual Green function given by 

= - i(T x0 A a (x°)4(x ro )) (21) 

= - t6(x Q - x' 5 )(A Q (x°)4(^°)) ± ^{x r ° - x°)(Al(x ,8 )A a (x°)}. 

The operators A a (x°) can be build out of any number of field operators (fermions and/or 
bosons) and their light cone time dependence in the Heisenberg picture follows Eq. (jSJ). 
For the free fermion field, i.e. A{x ) = ty(x°), the standard light-front propagator (jSJ) is 
recovered. 

^From the definition of the Green function at finite temperature we can derive the 
Dyson equation 

id-oG^-*' = 5(x 5 - x'«)N:l + R^'\ (22) 
where c\ is the derivation with respect to x° and 

= (K,4k(* 5 )>, <~*' 5 = -i(TAAa,H](x°)4(x ro )}. (23) 
With the following definition of the mass matrix 

M^ = J2[dx*lC 4 (G-^f (24) 

a' 

the Dyson equation (I2"2i) can be written as 

idoGtr' 5 = 5(x° - x' 5 )< + £ [ d^Mt^G%f . (25) 

The mass matrix describes the modification of the particles energy and of the interaction 
due to the medium. The expression (124"]) for the mass operator is not well suited for 
practical computations. However, after some formal manipulations, the mass operator 
reads 

= J2l^-^)(U a ,H],Al]± (x + )) 

a' 

- i( T x[ A a, H](x + )[H, A^,](x?)) irr .] (iV- 1 )^, . (26) 

This equation separates the mass matrix in a instantaneous term Mq^i related to the 
meanfield approximation and in a retardation (or memory) term M*~^ , i.e. 

M$T* = 5(x° - * 5 )M , a/3 , + Mf a f. (27) 

Here, we neglect the retarded part of the mass operator which leads to the following 
equation of motion 

id x o - Mf) G^~ x '° = 6 (x° - N x \ (28) 

where the indices a, (3 are suppressed for simplicity. 

For an explicit calculation of the expression for M and the normalisation factor iV 
assumptions about the medium are needed. Here, we assume homogeneous matter 
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of non-interacting quarks and antiquarks. In this quasi-particle approximation the 
calculation of the mass operator and of the normalisation factor is performed in 
momentum space by introducing creations operators for the particles s) and 
antiparticles d'(k,s). The assumptions imply directly the Fermi distribution functions 
f ± (k). For the one-body case one obtains the in- medium propagator ( 1551) from the Dyson 
expansion. To investigate mesonic bound states we have to calculate the normalisation 
factor of the quark- ant iquark system using (T25|) . It reads 



N 2 (k l ,k 2 ) = l-f+(k l )-fi(k 2 ) (29) 

and corresponds to the well-known Pauli blocking factor. Besides mesons three-quark 
bound states like the nucleon are subject of this paper, where we neglect the antiparticles 
degree of freedom. This is justified in light-front dynamics because pair creation 
processes are likely to be suppressed. The normalisation factor of the two-quark system 
is given by 

N 2 (k 1 ,k 2 ) = I-/1-/2, (30) 

where fi = / 4 + (fcj) for simplicity. The corresponding factor of the three-quark system is 

N 3 (h, k 2 , k 3 ) = h] 2 h - /1/2/3, (31) 

where we used the notation fi = 1 — fi- In general, the mass operator M can be 
separated into two terms: the first term, Hq, contains the self energy corrections in the 
effective mass and the second term describes the in-medium modified interaction. In 
the t-matrix equation we neglect the self energy term and consider in the following the 
term Hq only. 



The resolvent operator can be defined by 



< ] iz) = —L=r, z i a (Ho) , (32) 
z — n 

where n = 2, 3 indicates the two- and three-quark correlations respectively. Hence the 
Dyson equation f[2"51) yields 

(z) = irf (z) N n , (33) 

In this way the dominant medium effects due to Pauli blocking and self energy 
corrections are systematically included in the relativistic equations for the bound states, 
because the t-matrix equations contain G Q n ^ explicitly. In the case n = 2 the resulting 
equation for the two-body t-matrix has the same formal structure as the Feynman- 
Galitzkii equation 

T 2 = K + KRf ] N 2 T 2 (34) 

where K represents the two-body interaction kernel. For n = 3 the starting point is 
also given by the Bethe-Salpeter equation 

T 3 = K + KR^ ] N 3 T 3 . (35) 



Light-front NJL model at finite temperature and density 



8 



A derivation of relativistic three-body equation on the light-front is formally identical 
to the nonrelativistic [42, 43J case, if we neglecting antiparticle degrees of freedom. 
Following [42], 03j HI] the resulting equation for the vertex function \T a ) reads 

3 

|r a > = £(1 - S^N^R^Tp), (36) 

(3=1 

where N% , and R% a indicate the Pauli blocking factor, two-body t-matrix and the 
two-body resolvent in the channel /3. 



4. NJL model at finite temperature and density on the light-front 

In this Section we apply the light-front finite temperature field theory to the two flavour 
NJL model 

C = i)(i ltl d» - m )ij + G ((V^) 2 + (^75T</0 2 ) , (37) 

where r are the Pauli matrices. In the meanfield approximation the chiral quark 
condensate contributes to the dynamical generated mass m = mo — 2g{^jip). Using 
f l20|) to compute the chiral quark condensate leads to the in-medium gap equation 

/dk + d 2 k< 
fe+(2?r) > (l - f + (k + , k ± ) - f-(k + , k ± )). (38) 

The medium independent first summand in the integral in (|38|) needs to be regularised, 
as usual the temperature modifications to the mass m(T, /i) are finite. We like to 
use a consistent regularisation scheme for the gap equation and the related two-body 
calculations. Therefore the two-body bound state problem shall be discussed before 
treating the issue of regularisation of (l38l) . The dependence of the constituent quark 
mass on the medium has been calculated within the NJL model on the light-front in 
Ref. [H]. 



4-1. Pions 

The pion t- matrix T n (k) is obtained by the equation fl34|) that represented in momentum 
space reads 

T n (k) = K + J ^KG(q + k/2)G(q - k/2)T n {k), (39) 

where K is some irreducible interaction kernel and G(q) the quark propagator taking 
into account the thermal mass (1381) . To arrive at fl39|) we replaced the two-body Green 
function G ( o\q) by G{q - k/2)G(q + k/2) consistent with (1331). 

The separable interaction kernel in the pseudo-scalar channel reads 

K aPa s = -2iG(j 5 Ti) a p(jBTi) 7 s- (40) 
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Introducing a reduced t-matrix t n (k) via T n (k) a pys — (75 r j)o / 9^7r(^)(75 T 'j)7<5 the solution 
to ([39]) is 

M *> = T+mZm' (41) 

where U n (k 2 ) denotes the contribution of the bubble diagram for the pion. Evaluation 
of Il 7r (A; 2 ) using (1331) (or equivalently (120]) ) leads to 

dg^L M 2 2 |o(^ gO (1 ~ / + (M 2 |q) - /- (M 2 | )) 
x(1-x)(2tt) 3 M 2 Q (x,q ± )-k 2 l{ } 

which depends on the mass of the virtual quark- ant iquark system M^, (x,q±) = 
(q\ + m 2 )/x(l — x). We introduced the longitudinal momentum fraction x = q + /k + . 
The momentum dependence of the Fermi blocking factors can also be written in terms 
of M 2 |o as 



2|0J 



exp<(/3( } + l 



(43) 



The thermal mass of the pion m n (T,fi) is determined by the pole of t w (k). 

The integrals involved in (|38j) and (|42j) are divergent and need regularisation. We 
utilise the invariant Lepage-Brodsky (LB) cut-off scheme, which restricts the mass M 2 \o 
by 



M! l0 (x,q ± ) = ^^- ) <Al B . (44) 

We are left with three model parameters, namely the quark current mass mo, the 
coupling constant G, and the cut-off Alb- These are adjusted such that the pion mass 
m n = 140 MeV, the pion decay constant f n = 93 MeV, and constituent quark mass 
m = 336 MeV are obtained correctly One finds the values 

G = 5.51 eV, (45) 
m =5.67 MeV, 
A LB = 1428 MeV. 



These correspond to the case I discussed in [39] • Following Ref. [45] the gap equation 
with LB regularisation in light-front form and the gap equation with 3-momentum (3M) 
cut-off (k 2 < A| M ) in instant form are equivalent if one chooses a medium dependent 
LB cut-off 

A LB (T, /U ) = 4(A^ / + m 2 (T,/i)). (46) 

The medium dependence of the cut-off may seem artificial at first, but one should notice 
that the relevant energy scale of the effective field theory changes with the medium 
parameters. 

The in-medium pion mass for different values of the chemical potential is shown 
in Figure l.A. The Mott dissociation line is given by the intersection of m,r(T, fi) and 
the continuum 2m (T, /i). In the NJL model discussed the chiral phase transition occurs 




Figure 1. (A) The pion mass as a function of T for different /i. (B) The diquark 
mass as a function of T for different fi. The continuum is given by 2m. The lines of 
the two-body masses end at the Mott dissociation points. 



when m(T c , /i c ) = m(0, 0)/2 holds [39]. In Figure 3.B the chiral phase transition (T c , /i c ) 
is plotted as solid line and one reads off the critical temperature for vanishing chemical 
potential as T c = 190 MeV which is compatible with recent lattice calculations with 
almost physical quark masses [8]. The dashed line represents the pion dissociation. It 
worth noting that the pion dissociation and the chiral phase transition line are located 
in a narrow band of the phase diagram, which is of course expected by Goldstone boson 
character of the pion. However, both lines do not coincide in this simple model of quark 
matter, which means that deconfinement and chiral restoration happen at different 
temperatures and densities in accordance with certain lattice computations [46J. 



4-2. Diquarks and Colour Superconductivity 

The quark-quark interaction is constructed by Fierz transformation of fl3?j) . see 
Reference [47] for details. Let us consider scalar, isospin singulet, colour- antitriplet 
diquarks, i.e. we choose the following kernel 

K a ^ s = 2tG s ( 1 ,CT 2 X a ) a(S (C- 1 l5 T 2 X% 5j (47) 

with a = 2, 5, 7 for the anti-symmetric Gell-Mann matrices in standard representation 
and the charge conjugation matrix C = ij 2 ^ . Analogous to the pion case the reduced 
t-matrix td(k) is introduced and the solution 

2iG 
l + 2GJl d (k 2 ) 
contains the diquark loop 

n (h 2, _ f dx*q ± M 2 yx,gj)(l-2r(M 2|0 )) 
d[ ] ~ J x(1-x)(2tt)3 kqjx,q ± )-k* ' l49j 

LB 1 

where the normalisation factor (I3"0"j) is present. 
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Results for the medium dependence of the diquark mass m d are presented in Figure 
l.B. The coupling in the diquark interaction channel is fixed to G&jG^ = l/v2 which 
leads to the mass m d ~ 554 MeV of an isolated diquark. Interestingly the dependence 
of the in-medium diquark mass on temperature changes drastically when the chemical 
potential is increased. While for small /z the diquark mass decreases with temperature, 
the behavior is opposite for /i > 200 MeV. For even larger chemical potentials the 
chiral phase transition occurs and the diquarks cease to exist as bound states. This 
turn-around is caused by the competition of two effects, namely Pauli blocking and the 
in-medium quark mass. The contribution to the diquark loop free from Pauli blockings 
is still medium dependent and leads to a decreasing bound state mass as the quark mass 
reduces. Due to the different couplings in the pion and diquark channel one observes 
this differing behavior of the bound state mass even when the loop integrals ( 1421) and 
( 1491 are equal, e.g. for /i = 0. 

The interaction in the above discussed quark-quark channel is attractive. Therefore 
one expects the formation of a colour superconducting phase for sufficiently high 
densities and low temperatures. For two flavours the colour superconducting phase 
(2SC) is characterised by the condensate $ = (ip T C^T2\2ip) , cf. [40] for a review. The 
SU C (3) colour symmetry is broken down to SU(2). Usually one derives a gap equation 
for $ but here we follow a different strategy using the diquark t-matrix. The Thouless 
criterion 

m d (T, /i) = 2/i (50) 

enables us to compute the boundary to the colour superconducting phase jUj. Equation 
( |50l) is a condition for the condensation of bosonic diquarks but remains valid beyond 
the two quark threshold and leads to the cancelation of the pole present in (149]) for 
m d > 2m. Inserting (1501) into (]48l) the pole condition becomes 

1 _ f dx f d 2 q L M 2 2 |0 (x, q ± ) (l - 2/+(M 2 | )) 
2G d J x(l-x)J (2vr)3 M^x, q ± ) - Atf ' { } 

LB LB 

Our result is plotted in quark matter phase diagram Figure 3.B. The 2SC phase is located 
at temperatures below T < 92 MeV and for chemical potentials 300 MeV < \i < 702 
MeV. At T = 85 MeV and /i = 297 MeV the boundary of the superconducting phase and 
the chiral phase transition meet. These results are in general agreement with specific 
scenarios of more sophisticated instant form calculations including various condensates 
and colour superconducting phases [40] . 

4-3. Nucleon 

On the light-front the treatment of the spin is technically involved and for the time 
being we average over the spin projections. This procedure, used before in Ref. [49] to 
describe the proton electric form factor with a zero-range interaction, can be justified 
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in quark matter, because the spins can be regarded as washed out by the medium. 
Equation fl36|) evaluated in momentum space for a zero-range interaction reads 

n<l) = ^!M I d 4 kR?\k)N 2 (k,P 3 -q-k)4°\P 3 -q-k)r(k), (52) 



(2tt)< 

where are n-body free boson propagators, P^ is the three-body energy momentum 
vector in the center of mass system and the energy momentum vector of the two-body 
subsystem is given by P 2 = P£ — g M . The two-body scattering amplitude t(P 2 ) reads 

t(P 2 ) = (iX- 1 - B(P 2 ))~ l , (53) 

where the loop integral is given by 

i f dxd 2 k ± l-2f + (x,k 2 ± ) 
B[P2) ~ "2(2^ J x{l-x) P 2 2 -M 2 ■ (54) 



One notes the similarity of (H9l) and (1541) but the former integral accounts for the spins 
of the quarks which then lead to the M^ (x,q±) factor in the numerator. After the 
integration of fl52l) the integral equation for the vertex function can be written as 

T ^-W^Cw^) (55) 

where the arguments of the blocking factors read 

F(x, y; k±, q±) = f(x, k\) +f(l- x - y , (k ± + q ± ) 2 ). (56) 

The integral involving -6(^2) has a logarithmic divergence. In order to investigate a 
solution of the three-body bound state equation, we now extend the LB-regularisation 
explained in Section ETT1 to the mass of the virtual three-particle state M 3 | , which is the 
sum of the on-shell minus-components of the three particles, i.e. 

M 2 2 |0 , M 3 2 |0 < A 2 . (57) 

It is not meaningful to use the same value of the cut-off used in the two-body calculations, 
because the mass of the virtual three-particle state should be M 3 | > 3m and in this 
case we should have a too small integration range. Therefore, we express the parameter 
in units of the quark mass following [5UI EI], i-e. A = vm. The regularised two-body 
propagator is then modified as follows 



dxd 2 k ± l-2f+(x,k 



2 



5a(P2) ""W J x(l-x) P|-M? |0 {5 ' 



and the three-quark equation becomes 



rA(y^ = ^wf^ (59) 

a k± —2 —2 1 Ala;, k±). 



Mf |0 <A2 
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Figure 2. Binding energy of the three-quark bound state at various temperatures as 
indicated for A = 4m (left panel) and for A = 10 15 m (right panel) 



We have investigated the isolated two- and three-body bound states as a function 
of the coupling A for different values of the cut-off parameter and the stability of the 
three-body problem elsewhere [501 ED- Choosing m = 336 MeV for the quark mass we 
constrain our results to the nucleon mass. This determines the function A(A) and the 
model is parameterised by the cut-off A only. In the following we use the two limiting 
cases A = 4m and A = 10 15 m. 

The solution of allows us to calculate the three-quark binding energy at finite 
temperatures and the chemical potentials for the different cut-offs defined as 

B 3 {T, fi) = m(T, fi) + M 2B {T, //) - M 3B {T, //). (60) 

In Figure 2. the binding energy B3 as function of the chemical potential for constant 
values of the temperature using A = 4m and A = 10 15 m is shown. B% becomes smaller by 
increasing chemical potential for a constant value of the temperature. We can calculate 
at which temperature and chemical potential the binding energy vanishes and therefore 
the three-quark bound states disappear. As before the values of T and /x, for which 
the nucleon dissociation occurs, define the Mott lines and are shown on the left panel 
of Figure 3. in comparison with the chiral phase transition. At low temperatures the 
dependence on the cut-off is mild, but at zero density the different cut-offs A give 
significantly different values for the transition temperature. This can be understood by 
considering the important role of the colour screening at high density which cancels the 
details of the dynamics given by the different values of the cut-off. At low densities the 
effect of the colour screening is negligible and the different dynamics are perceivable. 
The Mott lines qualitatively follow the chiral phase transition given by the solid line 
and the nucleon dissociation occurs in the chiral broken phase. An overview of the 
phase diagram is given on the right side of Figure 3. The interval between these two 
transitions may be regarded as the confinement-deconfinement region. 
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Figure 3. (left panel) The Mott lines of the nucleon for A = 4m and A = 10 15 m. 
(right panel) The quark matter phase diagram in the light cone NJL model. The solid 
line shows the chiral phase transition. The dashed line is the Mott dissociation line of 
the pion and the dashed-dotted line is the Mott dissociation of the nucleon for A = 4m. 
Finally the transition between the deconfined and the 2SC phase is given by the dotted 
line. 



5. Conclusions 

Starting from a quasi-particle concept of quarks embedded in a hot and dense medium, 
we presented a unified description of the most important quark correlations present 
at the phase transition of QCD. Two challenges have been faced, correlations in a 
many-particle system and relativity. Correlations within a quasi-particle picture lead, 
in the simplest case, to (properly modified) few-body equations for two- and three-body 
states. Relativity leads to kinematical differences and the possibility of pair creation. 
In addition the strength of the interaction leads to a strong modification of the self 
energy of the particles. These challenges have been tackled utilising the light-front form 
of relativity. 

The basically Hamiltonian concept of the light-front form allows us to carry 
the Hamiltonian formulation of nonequilibrium quantum statistics along with the 
Dyson expansion (or cluster expansion) of many-body Greens functions through to the 
relativistic regime. For the sake of a unified description we have utilised a simple, but 
rather useful and hence widely investigated model of QCD, the Nambu-Jona-Lasino 
model. The basic ingredients of QCD relevant for finite temperature studies are included 
except for the confinement property. It is reasonable to assume confinement to be 
screened strong enough due to the presence of the other charged sources in the plasma. 
In fact, the major quality of this interaction is the zero-range property. Within this 
context we determine the in-medium properties of pions, diquarks and nucleons and 
investigated the whole temperature-density plane of the quark matter phase diagram 
including the chiral phase transition, the onset of the colour superconducting phase 
and the dissociation lines of pions and nucleons. The later ones give a measure were 



Light-front NJL model at finite temperature and density 



15 



the residual interaction between the quarks is weak enough to allow the onset of the 
confinement-deconfinement phase transition (assuming the colour dependent confining 
mechanism already screened accordingly). 

The physics field that can now be tackled is quite broad ranging from such intriguing 
questions how hadrons form during the early plasma phase of the universe, whether a 
colour superconducting phase exists, up to what the nature of QCD phase transition 
is. In addition, going beyond the quasi-particle approach seems to be straightforward 
in the Dyson expansion, however, it is still intriguing to arrive at self consistent 
results. Along this line a next step is the implementation of the full three-body it- 
matrix in a computation of the critical temperature of condensation/pairing of diquarks. 
The present formalism is well suited to address this very interesting problem, hardly 
investigated in literature. 

Besides the Dyson expansion in an equilibrated quasi-particle gas the light-front 
formulation allows us also to derive relativistic transport equations. First calculations 
that go beyond the NJL model have already been performed in 1+1 dimensions. Further 
effort is needed to approach the full 3+1 domain utilising e.g. transverse lattice methods 
or alike. 
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